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Abstract. In this paper we investigate tlie relationship between the mixing times of the Glauber dy- 
namics of a statistical mechanical system with its thermodynamic equilibrium structure. For this we 
consider the mean-field Blume-Capel model, one of the simplest statistical mechanical models that ex- 
hibits the following intricate phase transition structure: within a two dimensional parameter space there 
exists a curve at which the model undergoes a second-order, continuous phase transition, a curve where 
the model undergoes a first-order, discontinuous phase transition, and a tricritical point which separates 
the two curves. Wc determine the interface between the regions of slow and rapid mixing. In order 
to completely determine the region of rapid mixing, we employ a novel extension of the path coupling 
method, successfully proving rapid mixing even in the absence of contraction between neighboring states. 



1. Introduction 

The notion of mixing times is fundamental in the study of stochastic processes and their apphcations. 
In addition to its purely theoretical value, mixing times are widely used in computational and physical 
sciences. In computer science, mixing times represent running times for randomized algorithms such as 
Metropolis-Hastings that are responsible for providing high acctiracy solutions to a large class of problems 
where classical algorithms have impractical running times. In physics, one important use of mixing times 
is to determine the feasibility of simulating statistical mechanical systems via so called Glauber dynamics. 

An important question of mixing times for statistical mechanical models is its relationship with the 

thermodynamic phase transition structure of the system. More specifically, as a system undergoes an 
equilibrium phase transition with respect to some parameter; e.g. temperature, how does the correspond- 
ing mixing times behave? This question has been well studied with recent rigorous results for one class 
of statistical mechanical models, namely those that undergo second-order, continuous phase transitions 
like the famous Ising model [16] [15] [9]. For these models, it has been shown that the mixing times 
undergo a transition at precisely the thermodynamic phase transition point. However, for models that 
exhibit the other type of phase transition: first-order, discontinuous; e.g. Potts model with g > 2 [17], [8] 
and the Blume-Capel model [1], [2], [5]- [7], [14] with weak interaction, as far as well know, this question 
has not been previously addressed. In this paper, we prove that the mixing time transition does not 
coincide with the thermodynamic equilibrium phase transition for a model that exhibits a first-order, 
discontinuous phase transition. 

First-order, discontinuous phase transitions are more intricate than their counterparts, which makes 
rigorous analysis of these models traditionally more difficult. Furthermore, the more complex phase 
transition structure causes the models to fall outside the scope of standard mixing time techniques 
including the so-called 'path coupling' method [3]. The path coupling argument is a powerful tool used 
to prove rapid mixing for various Markov chains, including the Glauber dynamics for statistical mechanical 
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models that exhibit a second-order, continuous phase transition. However, due to the intricacy of a first- 
order, discontinuous phase transition, the standard path coupling argument can not be applied in all 
cases because contraction between couplings of all neighboring states does not exist. 

In this paper we prove the mixing time rates for the mean-field Blumc-Capel (BC) model [1], [2], [5]- 
[7], a statistical mechanical spin model ideally suited for the analysis of the relationship between the 
thermodynamic equilibrium behavior and mixing times due to its intricate phase transition structure. 
Specifically, the phase diagram of the BC model includes a curve at which the model undergoes a second- 
order, continuous phase transition, a curve where the model undergoes a first-order, discontinuous phase 
transition, and a tricritical point which separates the two curves. 

As mentioned above, in a subset of the first-order, discontinuous phase transition region of the BC model, 

the standard path coupling method does not apply as it docs throughout the second-order, continuous 
phase transition region. Therefore, we apply an innovative new method we call aggregate path coupling 
in order to prove rapid mixing for the BC model in the entire first-order, discontinuous phase transition 
region. While the standard path coupling method assumes contraction between every pair of states, our 
extended approach loosens this condition and requires contraction for only certain pairs of states. This 
new method is not restricted to only the BC model and can be applied to other statistical mechanical 
models that vmdergo a first-order, discontinuous phase transition and moreover is general enough to be 
a new tool in the theory of mixing times. 

The paper is organized as follows. The mean-field Blume-Capel model is introduced in Section 2. There 
we describe the equilibrium phase transition structure of the model and state the large deviation principle 
of the magnetization. In Section 3, wc define the Glauber dynamics for the BC model and the notion 
of mixing times of Markov chains. In Section 4, we introduce the path coupling for the BC model and 
prove the form of the mean coupling distance. In Section 5, we use the standard path coupling to prove 
rapid mixing in the single phase region where the BC model undergoes a second-order, continuous phase 
transition. In Section 6, we modify the path coupling argument to prove rapid mixing in the region where 
the BC model undergoes a first-order, discontinuous phase transition. The slow mixing results are stated 
and proved in Section 7. 



2. Mean-field Blume-Capel model 

Statistical mechanical models are defined in terms of the Hamiltonian function. For the mean-field 
Blume-Capel model, the Hamiltonian function on the configuration space Cl" = {—1, 0, 1}" is defined by 

" K ( V 

Hn,K{0j) = - - \ 

for configurations to = (wi, . . . Here K represents the interaction strength of the model. Then 

for inverse temperature /3, the mean-field Blume-Capel model is defined by the sequence of probability 

measures 

Pn,i3,Kiu}) = ^-^^--^exp[-/3if„,if(w)] 

where Zn{l3,K) = X^j^g^n exp[— is the normalizing constant called the partition function. 

In [14], using large deviations theory [11], the authors proved the phase transition structure of the BC 
model. The analysis of the Pn,/3,K was facilitated by expressing it in the form of a Curie- Weiss (mean-field 
Ising)-type model. This is done by absorbing the noninteracting component of the Hamiltonian into the 
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product measure P„ that assigns the probabiUty 3 " to each oj e il", obtaining 



(1) Pn,i3,K{d^) = - I ■ exp 



nl3K 



In this formula Sn{oj) equals the total spin Y^^=i^j^ Pn,p is the product measure on O" with identical 
one-dimensional marginals 

(2) Pp{dujj) = Yij^ ■ cxp(-/3w|) p{dujj), 

Z{(3) is the normahzing constant Jj^exp{—l3u)j)p{dcOj) = 1 + 2e~^, and Zn{l3,K) is the normalizing 

constant [Z(/3)]"/Z„(/3, ii'). 

Although Pn,i3,K has the form of a Curie- Weiss (mean-field Ising) model when rewritten as in (1), it 
is much more complicated because of the ^-dependent product measure P„^^ and the presence of the 
parameter K. These complications introduce new features to the BC model described above that are not 

present in the Curie- Weiss model [10]. 

The starting point of the analysis of the phase-transition structure of the BC model is the large deviation 
principle (LDP) satisfied by the spin per site or magnetization Sn/n with respect to Pn,0,K- In order to 
state the form of the rate function, we introduce the cumulant generating function of the measure 
defined in (2); for t S M this function is defined by 



c^{t) = log / exp(iwi)p^(rfwi) = log 

J A 



1 + 2e-^ 

We also introduce the Legendre-Fenchel transform of c^, which is defined for 2; e [—1, 1] by 

Jjsiz) = sup{tz - C/sit)} 

and is finite for z S [—1, 1]. Jjs is the rate function in Cramer's theorem, which is the LDP for Sn/n with 
respect to the product measures P„,^ [10, Thm. IL4.1] and is one of the components of the proof of the 
LDP for Sn/n with respect to the BC model Pn,fi,K- This LDP stated in the next theorem is proved in 

Theorem 3.3 in [14]. 

Theorem 2.1. For all P > and K > 0, with respect to Pn,p,K, Sn/n satisfies the LDP on [—1, 1] with 
exponential speed n and rate function 

Ip,K{z) = Mz) - /3Kz^ - mi{Jp{y) - pKy^}. 

3/eM 

In other words, for any closed subset F , 

(3) limsup - \ogPn,0,K{Sn/n & F] < -Ip,K{F) 
and for any open subset G, 

(4) liminf - \ogPn,ii,K{Sn/n e G] > -I^,k{G) 

where IfiMiA) = inf^g^ -?^£J,/f ( 



z 



The LDP in the above theorem implies that those z G [—1, 1] satisfying Ip^xiz) > have an exponentially 
small probability of being observed as n — )• 00. Hence we define the set of equilibrium macrostates by 

S0,K = {z€ [-1,1] ■I0,k{z) = O}- 
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For z e M we define 

(5) G0,K{z) = mz' -cp{2pKz) 

and as in [12] and [13] refer to it as the free energy functional of the model. The calculation of the zeroes 
of Ib.k cquivalcntly, the global minimum points of Jp^xiz) — ^Kz^ — is greatly facilitated by the 
following observations made in Proposition 3.4 in [14]: 

(1) The global minimum points of Jf).K{z) — PKz^ coincide with the global minimum points of G^,if, 
which are much easier to calculate. 

(2) The minimum values min^gRj J^^/f (z) — (iKz^} and xmn.z^^{Gp^K{z)} coincide. 

Item 1 gives the alternate characterization that 

(6) ^p,K = {z S [—1, 1] : z minimizes G^^^f (z)}. 

The free energy fimctional G^^x exhibits two distinct behaviors depending on whether P < Pc = log 4 or 
/3 > /3c- In the first case, the behavior is similar to the Curie- Weiss (mean-field Ising) model. Specifically, 
there exists a critical value Kc (13) defined in (7) such that for K < Kc (/3), G/s^k has a single minimum 
point at z = 0. At the critical value K = Kc Gjs^k develops symmetric non-zero minimum points 

and a local maximum point at z = 0. This behavior corresponds to a continuous, second-order phase 
transition and is illustrated in Figure 1. 




Figure 1. The free-energy functional G^,k for /3 < ft 

On the other hand, for > /3c, G^^k undergoes two transitions at the values denoted by Ki{^) and 
Kc^\p). For K < Ki{p), Gp^K again possesses a single minimum point at z = 0. At the first critical 
value Ki{P), Gp^K develops symmetric non-zero local minimum points in addition to the global minimum 
point at z = 0. These local minimum points are referred to as metastable states and we refer to Ki{P) 
as the metastable critical value. This value is defined implicitly in Lemma 3.9 of [14] as the unique value 
of K for which there exists a unique z > such that 

G"/3,Ki(/3)(2) = and G'f3,K^if3)i^) = 

As K increases from i^i(/3) to Kc^\/3), the local minimum points decrease until at K = Kc^\j3), the 
local minimum points reach zero and G^^k possesses three global minimum points. Therefore, for /3 > /3c, 

the BC model undergoes a phase transition at K = Ki^\l3), which is defined implicitly in [14]. Lastly, 

for K > Kc^\p), the symmetric non-zero minimum points drop below zero and thus Gpji has two 
symmetric non-zero global minimum points. This behavior corresponds to a discontinuous, first-order 
phase transition and is illustrated in Figure 2. 



MIXING TIMES FOR THE MEAN-FIELD BLUME-CAPEL MODEL VIA AGGREGATE PATH COUPLING 



5 




K< K,(fi) 



K = K,(|3) 



K,(P) < K < K'^'(P) 





Figure 2. The free-energy functional Gp^K for 13 > 



In the next two theorems, the structure of Sp^x corresponding to the behavior of G^ k just described is 
stated which depends on the relationship between /3 and the critical value /3c = log 4. We first describe 
^,3, if for < /3 < /3c and then for (5 > jic- In the first case £p^K undergoes a continuous bifurcation as K 
increases through the critical value Kc (/3) defined in (7); physically, this bifurcation corresponds to a 
second-order phase transition. The following theorem is proved in Theorem 3.6 in [14]. 

Theorem 2.2. For < P < (3c, we define 

if(2)(/3) = 



(7) 



4/3 



2^c^(0) 

For these values of /3, f has the following structure. 

(a) ForO<K< K^^\p), Ep,K = {0}. 

(b) For K > kP{(3), there exists z{l3,K) > such that £p,K = {±z{(3,K)}. 

f2l ('2') 

(c) z{l3,K) is a positive, increasing, continuous function for K > Kc (/3), and as 
z{f3^K) — > 0. Therefore, £p^K exhibits a continuous bifurcation at Kc 

(2) 

For /3 € (0,/3c), the curve {p,K'c'iP)) is the curve of second-order critical points. As we will see in a 
moment, for /3 G (/3c, oo) the BC model also has a curve of first-order critical points, which we denote by 

We now describe £p^K for /3 > /3c. In this case S/s^k undergoes a discontinuous bifurcation as K increases 
through an implicitly defined critical value. Physically, this bifurcation corresponds to a first-order phase 
transition. The following theorem is proved in Theorem 3.8 in [14]. 

Theorem 2.3. For all j3 > j3c, £p,K has the following structure in terms of the quantity Kc^\/3) defined 
implicitly for /3 > /3c on page 2231 of [14]. 

(a) ForO<K < K^^\p), E^^k = {0}. 

(b) There exists z{P, Kc^\p)) > such that £^ 



{0,±zi(3,Ki'\m 
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Figure 3. Equilibrium phase transition structure of the mean-field Blume-Capel model 



(c) For K > KP{f3) there exists z{f3, K) > such that S^^k = {±z{P, K)}. 

(d) z{j3,K) is a positive, increasing, continuous function for K > K^\j3), and as K ^ k''c'\I3)'^ , 
ziP,K) z{l3,K^^\/3)) > 0. Therefore, 

^I3,K exhibits a discontinuous bifurcation at Kc^\/3). 

The phase diagram of the BC model is depicted in Figure 3. The LDP stated in Theorem 2.1 imphes 
the following weak convergence result used in the proof of rapid mixing in the first-order, discontinuous 
phase transition region. It is part (a) of Theorem 6.5 in [14]. 

Theorem 2.4. For (3 and K for which 8p,K = {0}, 

Pn,i3,K{Sn/n € dx} => Sq as u ^ oo. 

We end this section with a final result that was not included in the original paper [14] but will be used 
in the proof of the slow mixing result for the BC model. The result states that not only do the global 
minimum point of Gp^K and Ip.K coincide, but so do the local minimum points. 

Lemma 2.5. A point z is a local minimum point of Gp^K if md only if it is a local minimum point of 

Proof. Assume that z is a local minimum point of Gp ^- Then z is a critical point of G^^k which implies 
that z ~ c'p{2/3Kz). By the theory of Legendre-Fenchel transforms, Jp{z) = (c^)^^(z) and thus 

IpMi?^) = JpiS) - 2f3KS - - 2/3KS = 0. 

Next, since 5 is a local minimum point of Gp^K, 

Gl^i^) > ^ 

Therefore, 

IIk{^) = 4{z)-2PK 



c"p{2l3Kz) < 



2I3K 



c'L{2f3K~z) 



2I3K > 
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and we conclude that z is a local minimum point of IjS^K- The other direction is obtained by reversing 
the argument. □ 



3. Glauber Dynamics and Mixing Times 



The Markov chain for statistical mechanical models studied in this paper is called the Glauber dy- 
namics, one of the most commonly used physical dynamics for these models. See [4] for more on Glauber 
dynamics. For the BC model, the Glauber dynamics evolve by selecting a vertex i at random and updat- 
ing the spin at i according to the distribution Pn,i3,Ki conditioned to agree with the spins at all vertices 
not equal to i. If the current configuration is u and vertex i is selected, then the chance of the spin at i 
is updated to +1 is equal to 

p2l}KS(u,i)/n 

(8) p+x{uj,i) :-- 



Q2l3KS{i^,i)/n _|_ gl3-(PK)/n _|_ g~2l3K S{uj ,i) /n 



where S{oj, i) = Ylij^ti is the total spin of the neighboring vertices of i. Similarly, the probabilities of i 
updating to and —1 are 



^l3-{PK)/n 



and 

^-2l3KSiuiA)/n 

^^^^ — ^2l3KS{uiA)/n ^ g/3-(^iK)/n ^ p-2l3K S{loA) /n 

p^i{ui, i) is increasing with respect to S{oj, i), i) is decreasing with respect to S{u!, i), and po{ui, i) 

is decreasing for S{uj,i) > and increasing for S{ui,i) < 0. 

The mixing time is a measure of the convergence rate of a Markov chain to its stationary distribution 
and is defined in terms of the total variation distance between two distributions jj. and u defined by 

\\H - i^Wtv = sup \n{A) - u{A)\ = ^ V \ii{x) - i/{x)\ 
Given the convergence of the Markov chain, we define the maximal distance to stationary to be 



d{t) = m&x\\P\x,-) - ttWtv 

where P*(a;, •) is the transition probability of the Markov chain starting in configuration x and tt is its 
stationary distribution. Then, given £ > 0, the mixing time of the Markov chain is defined by 

imix{£) = min{t : d{t) < s} 

In general, rigorous analysis of mixing times is very difficult and the proof of exact mixing time asymp- 
totics of even some basic chains remain elusive. See [16] for a detailed survey on the theory of mixing 
times. 



Rates of mixing times are generally categorized into two groups: rapid mixing which implies the mixing 
time exhibits polynomial growth with respect to the system size, and slow mixing which implies that 
the mixing time grows exponentially with the system size. Determining where a model undergoes rapid 
mixing is of major importance, as it is in this region that the application of the dynamics is physically 
feasible. 
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A classical tool in proving rapid mixing for Markov chains defined on graphs, including the Glauber 
dynamics of statistical mechanical models, is the path coupling technique [3]. It will be shown that this 
technique can be directly applied for the BC model in the second-order, continuous phase transition region 
but fails in a subset of the first-order, discontinuous phase transition region. For the latter region, we 
develop an extension of the path coupling method we call aggregate path coupling to prove rapid mixing. 
The path coupling method for the BC model is introduced in the next section. 



4. Path Coupling 

We begin by setting up the path coupling for the Glauber dynamics of the mean-field Blume-Capel 
model. Define the path metric p on O" = {—1, 0, 1}" by 

n 

(11) p(a,T)=^l{a, ^r,}, 

the number of sites at which the configurations a and r differ. 

Let a and r be two configurations with p{cf,t) = 1; i.e. a and r are neighboring configurations. The 
spins of fj and r agree everywhere except at a single vertex i. Assume that <Ji < Ti. We next describe the 
path coupling {X, Y) of one step of the Glauber dynamics starting in configuration a with one starting in 
configuration r. Pick a vertex k uniformly at random. We use a single random variable as the common 
source of noise to update both chains, so the two chains agree as often as possible. In particular, let U 
be a uniform random variable on [0, 1] and set 

)<U < p-i{(j,k) 

}-i{a,k) <U < p-i{a,k) +po{a,k) 




.i{T,k)+po{T,k) <U <1 



Since cr, < r^, for all j ^ i, S{a,j) < S{T,j) and thus 

p+i{T,k) > p+i{a,k) and p_i(r, fc) < p_i(<T, fc) 
The path metric p on the coupling above takes on the following possible values. 

f if k = i 

p{X, Y) = i 1 if k i and both chains updates the same 
[2 if k ^ i and the chains update differently 

The application of the path coupling technique to prove rapid mixing is dependent on whether the mean 
coupling distance with respect to the path metric p, denoted by Eo-,r [p{X, Y)] , contracts over all pairs of 
neighboring configurations. 



In the lemma below and following corollary, we derive a working form for the mean coupling distance. 
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Lemma 4.1. Let p be the path metric defined in (11) and {X,Y) be the path coupling where X and Y 
start in neighboring configurations a and r. Define 

2sinh(?^a;) 



(12) 
Then 



2cosh(?^a;) + e'5-^ 



Tl — 1 (tI — 1 ) 



Proof. Let n_i, no and n+i denote the number of —1, and +1 spins, respectively, not including the spin 
at vertex i, where the configurations differ, in configuration a. Note that n_i + no + n+i = n — 1. 

Define £(—1) to be the probability that X and Y update differently when the chosen vertex k ^ i is a 
— 1 spin. Similarly, define £(0) and £(+1). Then 



E,,r[p{X,Y)] 



^(1 - £(-1)) + ^(1 - £(0)) + ^(1 - £(+1)) 



+2 



£(-1) H £(0) + ^^£(+1) 

n n n 



n-1 n_i IX , , "+i r , t\ 

= 1 £(-1) H £(0) H £(+1) 

n n n n 

The probability that X and Y update differently when the chosen vertex fc 7^ Hs a — 1 spin is given by 
£(-1) = [P-i{cr,k) - P-i{T,k)] + [{p-i{cT,k) +po{a,k)) - {p-i{T,k) + po{T,k))] 
= \p+i{T,k)-p+i{a,k)] + \p-iia,k)-p-i{T,k)] 
= \p+iir,k)-p-i{T,k)] + \p-i{a,k)-p+i{a,k)] 

2sinh(g^(g„(r) + l)) 2sinh(g^(g„(a) + l)) 

2cosh(?^(5„(r) + 1)) + e^-^ 2cosh(?^(5„(a) + 1)) + e^" 

= ^0,K{{Sn{T) + 1)) - ^0MiSn{<T) + 1)) 

1 



Similarly, we have 
and 



V^,K{Sn{T)) - <Pl3,K{Sn{(T)) + O 

£(0) = V0,K{Sn{r)) - <ff>,K{Sn{cT)) 



£(+1) = ipp,K{{Sn{T) - 1)) - ^p,K{{Sn{'y) ' 1)) = (5„(t)) - ^0,K{Sn{<j)) + O 

Thus the proof is complete. 
For defined in (3), 



□ 



'Pli,K{x) = C, 

Thus we have the following corollary. 



a; (l + 0(l/n)) 



Corollary 4.2. Let p be the path metric defined in (11) and {X,Y) be the path coupling where X and Y 
start in neighboring configurations a and r. Then 

, , , - I - Cfl 1 2jjK^^ I +0 1 ^ 

n n ^ \ n 



E.,.[p(x,y)] = ^ + ("-^) 
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By the above corollary, we conclude that the mean coupling distance of a coupling starting in neighboring 
configurations contracts; i.e. E„^tIp{X,Y)] < p{cf,t) = 1, if 



2^K 



Snir) 



2PK 



2I3K 



213K 



< 



1 



1 



Since a and r are neighboring configurations and 5„(t) > Sn{a), this is equivalent to 



(13) 



2(5K 



Sn{cr) 
n 



< 



1 



2I3K 



Therefore, contraction of the mean coupling distance, and thus rapid mixing, depends on the concavity 
behavior of the function c^. This is also precisely what determines the type of thermodynamic equilibrium 
phase transition (continuous, second-order versus discontinuous, first-order) that is exhibited by the mean- 
field Blume-Capel model. We state the concavity behavior of in the next theorem which is proved in 
Theorem 3.5 in [14]. The results of the theorem are depicted in Figure 4 



Theorem 4.3. For > pc 

(14) 



log 4 define 



WciP) = cosh"^ {\^^ ~ ^^"^) - °" 



The following conclusions hold. 

(a) For < l3 < (3^, c'p{w) is strictly concave for w > 0. 

(b) For P > (3c, c'l^iw) is strictly convex forO < w < Wc{l3) andc'^{w) is strictly concave for w > Wc{(3). 

By part (a) of the above theorem, for /3 < Pc, c'^ix) < c^(0) = l/{2pKP {/3)). Therefore, by (13), the 
mean coupling distance contracts between all pairs of neighboring states whenever K < Kc {(3)- 



By contrast, for (3 > /3c, we will show that rapid mixing occurs whenever K < Ki{l3) where Ki{(3) is the 
mctastablc critical value introduced in Section 2 and depicted in Figure 2. However, since the supremum 
sup[_i 1] c|^(a;) > 2/3Ki{i3) ' having K < Kx{f3) is not sufficient for (13) to hold. That is, K < Ki{l3) 
does not imply the contraction of the mean coupling distance between all pairs of neighboring states. 
Regardless, we will prove rapid mixing for all K < K\{j3) in Section 6 using an extension to the path 
coupling method that we refer to as aggregate path coupling. 



We now prove the mixing times for the mean-field Blume-Capel model, which varies depending on the 
parameter values {j3, K) and their position with respect to the thermodynamic phase transition curves. 



5. Rapid Mixing for Continuous Phase Transition Region 



We begin by stating the standard path coupling argument used to prove rapid mixing for the mean- 
field Blume-Capel model in the continuous, second-order phase transition region. The result is proved in 
Corollary 14.7 of [16]. 

Theorem 5.1. Suppose the state space D, of a Markov chain is the vertex set of a graph with path metric 
p. Suppose that for each edge {a, r} there exists a coupling {X, Y) of the distributions P{a, •) and P{t, •) 
such that 

E„^r[p{X, ¥)] < p{a, r)e-" for some a>0 

Then 

log(e) -|- log(diam(ri)) 



imix (^) — 



a 
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Figure 4. Behavior of c'g{w) for large and small /3. 



In this section, we assume P < Pc which imphes that the BC model undergoes a continuous, second-order 
phase transition &t K = Kc (/3) defined in (7). By Theorem 4.3, for P < Pc, c'p{x) is concave for x > 0. 
See the first graph of Figure 4 as reference. We next state and prove the rapid mixing result for the 
mean-field Blume-Capel model in the second-order, continuous phase transition regime. 

Theorem 5.2. Let tmixi^) be the mixing time for the Glauber dynamics of the mean-field Blume-Capel 

(2) 

model on n vertices and Kc (P) be the continuous phase transition curve defined in (7). Then for 
/? < ;9c = log4 and K < kP{P), 

n 



a 



-(logn + log(l/e)) 



for any a G ^0, ~~]^(2r~) '^'^'^ large enough. 



Proof. Let {X, Y) be a coupling of the Glauber dynamics of the BC model that begin in neighboring 
configurations a and r with respect to the path metric p defined in (11). By Corollary 4.2 of Lemma 4.1, 



,1 (n-l) 



2PK^)-C'J2PK^ 



O 



Observe that is an even function and that for P < Pc, supc^(x) — c^(0). Therefore, by the mean value 

X 

theorem and Theorem 2.2, 

E,MX,n < • |l-("-l)<Wn)ci;(0)l 

< 



< 
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/ '—K\ (2) 

for any a S (0, ^(2) 1 and n sufficiently large. Thus, for K < Kc (/3), we can apply Theorem 5.1 
where the diameter of the configuration space of the BC model Cl" is n to complete the proof. □ 

6. Rapid Mixing for Discontinuous Phase Transition Region 

Here we consider the region /3 > /3c, where the mean- field Blume-Capel model undergoes a first-order 
discontinuous phase transition. In this region, the function c^(x) which determines whether the mean 
coupling distance contracts (Corollary 4.2) is no longer strictly concave for .t > (Theorem 4.3). See 
the second graph in Figure 4 for reference. We will show that rapid mixing occurs whenever K < Ki{p) 
where Ki{/3) is the metastable critical value defined in Section 2 and depicted in Figure 2. 

As shown in Section 4, in order to apply the standard path coupling technique of Theorem 5.1, we 
need the inequality (13) to hold for all values of S'„((t) and thus sup[_]^ ^ c'l^i^) < 2/51? • However since 
sup[_i 1] c'^{x) > 2/3 if 1 ' condition K < i^i(/3) is not sufficient for the contraction of the mean 
coupling distance between all pairs of neighboring states which is required to prove rapid mixing using 
the standard path coupling technique stated in Theorem 5.1. 

In order to prove rapid mixing in the region where (3 > (3c and K < Ki{l3), we take advantage of the 
result in Theorem 2.4 which states the weak convergence of the magnetization S„/n to a point-mass 
at the origin. Thus, in the coupling of the dynamics, the magnetization of the process that starts at 
equilibrium will stay mainly near the origin. As a result, for two starting configurations a and r, one 
of which has near-zero magnetization {Sn{a)/n w 0), the mean coupling distance of a coupling starting 
in these configurations will be the aggregate of the mean coupling distances between neighboring states 
along a minimal path connecting the two configurations. Although not all pairs of neighbors in the path 
will contract, we show that in the aggregate, contraction between the two configurations still holds. 

In the next lemma we prove contraction of the mean coupling distance in the aggregate and then the 
rapid mixing result for the mean-field Blume-Capel model is proved in the theorem following the lemma 
by applying the new aggregate path coupling method. 

Lemma 6.1. Let {X,Y) he a coupling of the Glauber dynamics of the BC model that begin in config- 
urations a and t, not necessarily neighbors with respect to the path metric p defined in (11). Suppose 

j3 > j3c and K < Ki{j3). Then for any a S ^0, ^j^f^^ )^ ^ there exists an e > such that, asymptotically 
as n ^ 00, 

E,,<,[p(X,y)]<e-«/Xa,r) 

whenever \Sn{(T)\ < en. 

Proof Observe that for (3 > l3c a.nd K < Ki{l3), 

We will show that for a given a' e (^ 20Ki{fi) ' 1^)' t^i^re exists £ > such that 

(15) c'l^i^) ~ c'l^i^o) < oi'{x — xq) whenever |a;o| < e 

as c.g{x) is a continuously diff'erentiable increasing odd function and c'JO) = 0. 
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In order to show (15), observe that c"p{Q) = ^^j^m < 2pKi {/3) ' ^^'^ since c'^ is continuous, there exists 
a ^ > such that 



c'^{x) < a' whenever \x\ < d 



The mean value theorem imphes that 



c^(a;) - c'/^ixo) < a'{x - xq) Vxo,a; e {-S,6) 



Now, let e = a'+i/(2^Xi(^j) '^ < ^- Then for any |a;o| < £ and \x\ > 6, 



\r' / IK 1^1 + l^ol < + _ k-^ol 1 + \x-xo\ l + e/6 

mx) C0[xo)\ _ 2pK,{l3) - 2pK^{p) 2pK,{p) ' |1 - xo/x\ " 2/3/^i(^) ' l-e/6 



a \x — xo\ 



Without loss of generality suppose that S'n(cr) < S'„(t). Let {a = xq^xi, . . . ,Xr = t) be a path connecting 
(7 to r and monotone increasing in p such that (xj-i, Xj) are neighboring configurations. Here r = p{a, r). 
Then by Corollary 4.2 of Lemma 4.1 and (15), we have for |5„(cr)| < en and asymptotically as n — )• oo, 



E,,,[p(x,y)] < ^E,._„,jp(Xi_i,xo] 

(n-1) , , (n-1) 



< 



n 

(n-1) 



, f2pK 



n 



Sn{r) - d 



, (2PK 



+ p((T,r)-0 



p a,T + ^ -{Sn{r) - Sn{a))^ + p{(t,t) -Oi^ 

n n n 



1 - 



1 - 2pKa' 



O 



1 



< e-"/>(a,T) 



This completes the proof. 



□ 



Theorem 6.2. Let t„iix{e) be the mixing time for the Glauber dynamics of the mean-field Blume-Capel 
model on n vertices and Ki{P) be the metastable critical point. Then, for j5 > j5c and K < Ki{j3), 



tmix{£) < -(logn + log(l/£)) 
a 



for any a € (^0, ^KiI/s)^ ^ '^'"'^ large enough. 
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Proof. Let {Xt,Yt) be a coupling such that Yq '''= Pn,i3,K, the stationary distribution. For a given 
^Kiifi)^ ) ' ^ ^® ^ Lemma 6.1. Then for sufficiently large n, 

\\P\Xo,-)-Pn,0,K\\TV < P{Xt^Yt} 

= P{p{Xt,Yt) > 1} 

< E[p{Xt,Yt)] 

< E[p{Xt,Yt) I \Sn{Yt-i)\ < sn] ■ P{|5„(Ft_i)| < en} + 2nP{|S„(yt_i)| > sn} 

< e-«/"E[p(Xt_i,y(_i);|5„(yt_i)| <£n]-P{|5„(yi_i)| <£n} + 2nP{|5„(Ft_i)| >£n} 

< e-«/"E[p(Xt_i,yt_i)] • P{\Sn{Yt-i)\ < en} + 2nP{\Sn{Yt-i)\ > sn} 

t-1 t-1 

< e-"*/^E[p{Xo,Yo)] • n^{l^n(n)l < en} + 2nJ2P{\Sn{Ys)\ > en} 

s=0 s=0 

< e-"*/"E[p(Xo, Yo)] ■ iPaj3,K{\Sn/n\ < £})* + 2ntP„j3^K{\Sn/n\ > e} 

< 2ne-"*/" • iPn^p,K{\Sn/n\ < e}f + 2ntP„,^,K{|5„/n| > e} 
We recall the result in Theorem 2.4 that for P > jic and K < Ki{j3) 

Pn.p.K{Sn/n G dx} 5o as n — > oo. 
Moreover, for any 7 > 1 and n sufficiently large, the LDP stated in Theorem 2.1 implies that 

\\P\Xo,-)-Pn,p,K\\TV < 2ne-'''/^ ■ {Pn,p,K{\Sn/n\ < e})* + 2ntPn,p,K{\Sn/n\> e} 

< ne""*/" • (1 - e-'^"^^'^(^)) ' + ine'^^'^'^^^^ 
This completes the proof. □ 



7. Slow Mixing 

We complete our mixing time analysis of the mean-field Blume-Caple model by determining the slow 

mixing region of the (/?, K) parameter space for both the continuous, second-order and discontinuous, 
first-order phase transition regions. To prove slow mixing for the BC model, we use the bottleneck or 
Cheeger constant argument which we state in the theorem below. The theorem is an application of 
Theorem 7.3 in [16] to the BC model. 

Theorem 7.1. For two configurations u and t, define the edge measure Q as follows: 
Q{co,t) = Pn,i3,K{io)P{u,T) and Q{A,B)= ^ Q(w,r) 

uj£A,t€B 

Here P(u!, r) is the transition probability corresponding to the Glauber dynamics of the mean-field Blume- 
Capel model. The bottleneck ratio of the set S is defined by 

Pn,fi,K{S) S:P„,p,k(S)<1 

Under these assumptions, we have 
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The next lemma states that bottlenecks, and thus slow mixing, occurs for the BC model whenever the 
function G^^k has a positive minimum point. 

Lemma 7.2. Suppose the free energy functional G^^k has a positive minimum (either local or global) 
point. Then there exists a positive constant b and a strictly positive function r{(3,K) such that 

t ■ > be''('3,A:)n 

Proof. Suppose G/j.a: has a minimum (cither local or global) point at 2 > 0. Furthermore, let z' be where 
Gp^K has a local maximum point such that < z' < z. Define the bottleneck set 

A = {uj: z' < Sn{oj)/n < 1} 

Since z' > 0, by the symmetry of Gp^x, Pn,i3,K{A) < \. With respect to the Glauber dynamics, in order 
to leave the set A in one step of the chain, the current configuration must lie in one of the two boundary 
sets 

Ai = {uj: =z' + -\ or A2 = {w. =z' + -\ 

nj y n n ) 

We take n sufficiently large so that z' + 2/n < z. 

For any 7 > 1, the large deviations upper bound (3) implies that the edge measure Q for the bottleneck 
set A satisfies 

Q{A,A') < Pn,p,K{Ai) + Pn,p,K{A2) < e'^^^'^^^' + S) 

for n large enough. Moreover, since the LDP implies that P„_^^ii-(^) > e""'''^'^'^^^^ the bottleneck ratio 
satisfies 

all for a given 7 > 1 and n large enough. Lastly, by Lemma 2.5 z is a local minimum point of Ip,K since 
it is a local minimum point of G/^^k, and thus Ii}^k{z') > Ii3^k{z) and the proof is complete by Theorem 
7.1. , , , ^ 

Next we state the slow mixing result for the mean-field Blume-Capel model. The result follows from 
Lemma 7.2 and the region in the (/3, K) parameter space at which the free-energy functional G^^k 
possesses a positive minimum point. These regions were determined in [14] and depicted in Figures 1 
and 2. 

Corollary 7.3. Let t^ix = tmixi^/^) be the mixing time for the Glauber dynamics of the mean-field 
Blume-Capel model on n vertices. For (a) P < Pc o-nd K > Kc and (b) P > Pc o-nd K > Ki{P), 
there exists a positive constant b and a strictly positive function r{P,K) such that 

t ■ > be'^^'^'-^'" 

^mix — 

We summarize the mixing time results for the mean-field Blume-Capel model and its relationship to the 

model's thermodynamic phase transition structure in Figure 5. As shown in the figure, in the second- 
order, continuous phase transition region (/3 < Pc) for the BC model, the mixing time transition coincides 
with the equilibrium phase transition. This is consistent with other models that exhibit this type of phase 
transition. On the other hand, this is not the case in the first-order, discontimious phase transition region 
(/? > Pc) where the mixing time transition occurs below the equilibrium phase transition. 
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Figure 5. Mixing times and equilibrium pliase transition structure of the mean- field Blume-Capel model 
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